Cumulative structure function in terms of 
nucleonic wave function of the nucleus 



M.A.Braun, Bv.M.Suslov land B.Vlahovic 
NCCU,Durham, NC, USA 



(N 
O 
O 
(N 

: Abstract 

m 

The structure function of the nucleus in the cumulative region x > 1 is studied in terms 
of nucleon degrees of freedom. At high the resulting expressions are presented as a sum 
of contributions from few- nucleon correlations. Two- nucleon correlations are studied in some 
. detail. Spin variables are averaged out. In the region 1 < x < 2 the structure functions 

are calculated for the relativistic interaction proposed by F. Gross et al. They are found to 



o 
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' fall with X faster than the exponential. For Carbon at x = 1.05, where the method is not 

^ : ngorou* applicable, they turn out to be rougly twice larger thau the experimental data, 

1 Introduction 

^ ' The cumulative phenomena, i.e. processes with nuclei in the kinematical region prohibited for 

C ■ the colliding free nucleons, are of considerable interest both from the practical and theoretical 

poits of view. On the practical side they give a possibility to effectively raise the energy of 
KJi \ the colliding particles, allowing for subthreshold production of high mass partciles. On the 

^ • theoretical side they allow to study small internucleon distances and thus study the nuclear 



matter at high densities and possibly its quark structure. One may speculate that in these 
phenomena the cold quark-gluon plasma is formed in the region of many overlapping nucleons. 

Most of the experimental data on the cumulative phenomena refer to particle production 
in the region outside the free nucleon kinematics. Obtained more than twenty years ago, 
they were actively discussed at that time. Many models were proposed to explain the data. 
However no general consensus was reached. Even the basic question, whether the quark 
rather than nucleon degrees of freedom were nesessary to describe the experimental results 
remained unsolved. The reason for this lies in the complexity of the production process. 
It inevitably involves soft interaction of the colliding particles, the description of which is 
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highly model dependent and thus introduces uncontrollable elements into the predictions for 
cumulative particle production rates. 

From the theoretical point of view a much cleaner observable is the structure function of 
the nucleus in the cumulative region, that is, at the values of x greater than unity. Unfortu- 
nately the relevant experimental information is quite scarce, limited to data on the deuteron 
[1] and carbon [2] at relatively low values of x < 1.5. However even these modest data present 
a clear possibility for theoretical analysis and testing of different approaches. A description 
of the cumulative structure function of the nucleus based on the QCD approach and quark 
structure of the nucleus was proposed in [3]. Potentially this approach can relate the ob- 
served magnitude and behaviour of the cumulative structure function to the fundamental 
QCD parameter A. However the results obtained in [3] strictly speaking refer to the extreme 
cumulative kinematical limits. They also result infrared unstable and require introduction of 
an infrared cutoff ("quark mass") as a parameter. With this single parameter both the slope 
and the magnitude of the cumulative structure function were described reasonably well in 
[3]. 

For all that, it remains unclear whether the description of the nuclear structure function 
at moderate cumulativity, that is for x ~ 2 3 necessarily requires using the quark language 
and thus throwing away most of the information accumulated from the study of the nucleus 
in terms of hadrons. It is worth mentioning that many years ago Frankfurt and Strikman 
strongly advocated the thesis that the cumulative phenomena could be well described in 
terms of nucleons, provided the relevant wave functions be appropriately relativized [4]. 

In this paper we pursue exactly this approach. Exploiting the fact that forms for the 
relativistic internucleon interaction have lately been proposed, which lead to excellent agree- 
ment with all low energy data, we use this interaction to calculate the cumulative nuclear 
structure function. It is important that the cumulative structure function is directly related 
to the internucleon interaction, unlike most of other properties which rather involve all sort 
of averages of the interaction. Thus our results constitute a stringent test on the proposed 
relativistic interaction. In principle our approach allows to find the cumulative structure 
function at all values of x < ^. However calculational difficulties grow very fast with x. For 
that reason in this paper, apart from presenting the general formalism, we calculate the nu- 
clear structure function only in the interval 1 < x < 2 from the contribution which is known 
as "two-nucleon correlation" [4], leaving the three-(and more-)nucleon correlations for future 
studies. 

The paper is organized as follows. In Sees. 2 and 3 we present our formalism and separate 
contributions from correlations of several nucleons. In Sec. 4 we analyze the two-nucleon 
correlation first for the simplified scalar case and then for the relaistic spinor case. Sec. 5 
presents our nucmerical results. Some conclusions are drawn in Sec. 6. In the Appendices 
the parameters of the used relativistic interaction are presented and also some calculational 
details. 
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2 Formalism 

At high values of the nuclear structure function F^^^ix^Q"^^ is given by the "impulse 
approximation" diagram shown in Fig. 1, corresponding to the leading twist contribution. 
Our notations are clear from Fig. 1. We define in the standard manner 

x = ^ 
2qp 

Note that x is defined respective to the collision with a nucleon. The mass of the nucleus is 
M = A{m - e), so that 

= — 2me = — 

In future we shall neglect the binding energies wherever it is possible. The physical region 
for the reaction is evidently 

{q + Apf > or 0<x<A 

All the region 

l<x<A (1) 

is cumulative. 

The separation of the structure function from the diagram shown in Fig.l is standard. 
Let VF/li^^ be imaginary part of the amplitude (the discontinuity in s divided by 2i). Then 
one presents 

w!.^' = {-9.^ + '-f) Pi^\-^Q') + ^ - ) - ) Pi^\-,Q') (2) 

with P = Ap the nucleus 4-momentum. Functions F^^^'^ are the nuclear structure functions. 

In the following we concentrate on ^2^^^ and suppress subindex 2. It is convenient to choose 
a system in which = 0. Then taking the ++ component of (2) one gets 



wi'i = ^F(^\x,Q') (3) 

This relation serves to read the structure function directly from the diagram shown in Fig. 
1. 

To move further for simplicity we assume nucleons to be spinless particles. However the 

final results will also be valid for realistic nucleons after averaging over spins. Some details 

(A) 

concerning inclusion of spins will be presented in Section 4. We present W^i^ as an integral 
over the momentum k of the active nucleon 

Wi^^ = J ^W^u\Q,k) 2Im^p,k) (4) 

where $ denotes the lower blob in Fig. 1 with the propagators of the active nucleon included 
and wliu'' is the imaginary part of the upper blob. The latter can be presented as in (2), 
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with the off-sheh structure functions F^^\x' , Q'^ , k'^) where x' = Q'^/{2qk). At high Q'^ at 
lower twist these structure functions do not depend on k"^ and coincide with the observable 
structure functions of the nucleon. At high in the system g+ = we have g_ ~ Q^, 
q_L ^ Q « q+ so that qp ~ q^p+ and qk ~ q-k^. So if the scaling variable of the active 
nucleon is z = then x' = xjz. With these relations we find from (3) and(4): 

F(^)(x,Q^) = \\yzzF^^) (^,Q2) / ^±|-^2Im$(p,fc) (5) 

where we used the fact that the structure function of the nucleon depends only on the scaling 
variable of the active nucleon. 

The integral over /c_ , kx. of the Im # can be related to the distribution in x of the nucleons 
in the target nucleus. To do this consider the baryonic form-factor B^^q) at q = described 
by the diagram shown in Fig. 2. In it the lower blob is the same forward scattering amplitude 
A;) as in Fig. 1 for the structure function, only non-cut. Blob $ is a function of two 
variables: the active nucleon virtuality rv? — k^ and the square of the cm. energy of the 
recoil particles s = {Ap — k)^. Neglecting the dependence of the nucleon form- factor on the 
nucleon virtuality we have: 

In terms of integration variables in the second integral 

e='Ml--^s--^kl (7) 
A ^ A - z A - z ^ ^ 

and the integration over s goes along the Fcynman path. 

In the complex s plane the singularities of $ come, first, from the standard unitarity 
righthand cut and, second, from the singularities in k"^. These lie along the positive values 
of k'^. However due to a coefficient in (7) they transform into a lefthand or righthand cut in 
s depending on whether the ratio z/(A — z) is positive or negative. If the cut coming from 
the singularities in k^ lies to the right, we shall get zero, since the integration contour in s 
then can be closed in the upper hemisphere with no singularities inside. So, as expected, a 
non-zero result only follows if 

0<z<A (8) 

when the cut from the singularities in fe^ lies to the left, and the result can be obtained by 
closing the contour around the unitarity cut. Taking the component of (6) we get 

7 — ^ — n r/ — Im$(s,A;^) = (9) 

(27r)2 Jq 2n{A - z) Js^ 2t: ^ ' ' ^ ' 

This relation can be rewritten as a normalization condition for the relativistic nuclear density 
in z and k± 

dz f d'^k± , , 9 , , , 
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where the density p{z, kjj is defined as an integral over the "-" component of k: 

Note that on the formal level, 

2lm^{s,k^) = j d\e-'^''{A\N\r)N{Q)\A) = 

J2(.27r)H\Pa + k - P){A\N^O)\A - 1, a){A - 1, a\N{0)\A) (12) 

a 

where \A) is the state of the nucleus at rest, |A— 1, a) are the recoiling nuclear states with A—1 
nucleons, including states with A — 1 free nucleons, and N{x) is the nucleon field operator. 
The amplitudes {A — l,a\N{0)\A) describe splitting of the nucleus into a virtual nucleon and 
recoil states. The summation over all states in (12) is restricted by the energy conservation 
law (the conservation of "-" component of the momentum in light-cone variables). Integration 
over k- lifts this restriction and, as is evident from (12), puts the operators N and at 
equal "time" = 0. This allows to demonstrate that in the light-cone formalism with a 
conserved number of nucleons p{x,k'j_) is indeed the relativistic nuclear density 



\ip{x,kr,X2,k2±; ...;xA,fcA,±)P (13) 

where ip is the standard light-cone wave function of the nucleus, symmetric in all the nucleons 
and normalized according to 



/ 



The derivation is presented in Appendix 1. 

The light-cone wave function ^ can be related to the non-relativistic wave function by 
taking the non-rclativistic limit in the normalization condition (14), which corresponds to 
Xj — > 1. In this limit (14) gives 



/ 



A 

\{(fkj\ip{x,kr,X2,k2r,,,,\XA,kA,i.)\^ = (167r3)V^-i (15) 

J=2 



where we introduced kjz = m{xj — 1) for j > 2. Prom this we conclude that in this limit 

V'(x, kr, X2,k2±; ,,,;xA, kA,±) ^ (167r3)^/W^-i)/V(A;i, fes, ^a) (16) 



where (f) is the non-relativistic wave function of the nucleus at rest. On the right-hand side 
kj denote the 3-dimensional momenta, their sum equal to zero. 
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Returning to the structure function we find from (5) 

F(^\x, Q') = aJ^^^ p{z) (17) 

where the distribution in z is just 

P(^) = / §^Pi^^kl) (18) 

normahzed according to 

/ fK^) = l (19) 
Jo z 

3 Cumulative vs. non-cumulative regions. Correlations 

As indicated, blob <I>(p, /c) includes two external propagators for the active nucleon and so 
has a double pole at k"^ = m?. If we neglect the binding, then with all recoil particle at rest 
we shall have k = p and so fc^ = m?, which implies that the distributions p{z, kj^) and p{z) 
are strongly peaked at z = 1 and k± = 0. These values correspond to the non-relativistic 
domain, where the density p is more or less known from the existing nuclear data at relatively 
small transferred momenta. 

Going back to the nuclear structure function we observe from Eq. (17) that one can 
distinguish between two different kinematical situations. If x < 1 (non-cumulative region) 
then the integration region includes the vicinity of z = 1, where as we have just seen, p[z) 
has a sharp maximum. Then, as a first approximation, one can take F{x',Q'^) out of the 
integral at this point (i.e. at x' = x) and integrate the rest using the normalization condition 
(19), which gives 

f(^Hx,Q^)=AF<^^\x,Q^) (20) 

One can improve this result by taking some form for p(z) and integrating over z in a proper 
way. Then one will get corrections to (20) (the EMC effect). It is worth noting that it is always 
possible to match the experimental data by appropriately choosing p at finite (relativistic) 

values of I z — 1 1 . 

However in the cumulative region x > 1 the point z = 1 stays outside the integration 
region in (17). This means that in this interval of x values of p are important which come from 
$(p, k) with the active nucleon far off-shell, well outside the non-relativistic region. Prom the 
point of view of the latter, it is the far relativistic asymptotics of the nuclear dynamics which 
now matters. 

The technique to study this asymptotics is in fact standard. Let us consider a simple 
example of the deuteron, for which <^{p.k) reduces to the square of the Bethe-Salpiter wave 
function 'ijj^'^\k, 2p — k) satisfying the equation 

{n?-ki){m'-kl)^^<'\k^M) = I ^Vik,,k2;k[,k'^)i,(^\k[,k'^) (21) 
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where ki + k2 = k'l + k2 = 2p. In the vicinity of the mass shell kf = k2 = this 
equation reduces to the ordinary Schroedinger equation and correspondingly ip^'^^ reduces 
to its nonrelativistic limit V'/vR- Now suppose wc want to study the asymptotics of ip^"^^ far 
off shell. A simple way do to it is to rewrite (22) as 

^^'Hki,k2) = , 2 .2I 2 .2, / Pkvik^,k2■,k[,k'2)^P^''\k[,k'2) (22) 
(m"^ — ki){m'^ — kI) J (27r)^ 

and take into account that the integration is dominated by the close to on-shell values of the 

momenta, that is, by the non-relativistic region. So inside the integral one can substitute 
V'^'') by its non-relativistic limit and moreover take the interaction V outside the integral at 
small values of the 3-momenta k'12 corresponding to the nonrelativistic region. In this way 
one gets the desired asymptotics as 

~ (m^ - - fca '^C-'-fe^"'^)/ (23) 

where p = (m — e,0). In terms of diagrams this procedure corresponds to transforming the 
NND vertex as shown in Fig. 3, where it is assumed that the internal lines correspond to the 
non-relativistic deuteron and the external ones correspond to the relativistic deuteron. 

For a nucleus consisting of more than two nucleons this procedure leads to the diagrams 
shown in Fig. 4a. If more than two nucleons are relativistic, it can be repeated, leading 
to diagrams like Fig. 46. Following [4] one speaks of two- (Fig. 4a) or three- (Fig. 46) or 
more-body correlations in the nucleus to describe two, three or more relativistic nucleons. 

From the pure kinematical considerations one immediately concludes that the A;-fold cor- 
relation gives a contribution to the distribution p{x) in the region 

< X < A; 

Equivalently, choosing a particular interval of x 

n-l<x<n (24) 

the contribution to the nuclear structure function in this region comes from correlations of k 
nucleons with k > n. 

At large virtualities, with interaction and propagators falling like powers, one finds that 
the contribution from the /c-fold correlations diminishes fast with the growth of k. The bulk 
of the contribution in the region (24) will then come from precisely the n-fold correlation, the 
ones with k > n giving only a small correction. So, in the first approximation, to study the 
structure function in the region (24) one has to take into account only the ra-fold correlation. 

Up to now the cumulative nuclear structure function has been studied experimentally 
only in the region 

1< X < 2 (25) 

This gives motivation to start from the pair correlations in the nucleus, which will be the 
subject of the next sections. 
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4 Two-body correlations 
4.1 Scaletr nucleons 

To clearly formulate our approach we continue with the scalar nucleons case. As explained 
in the previous section, the contribution of the 2-body correlations to function cor- 
responds to the diagram shown in Fig. 5 in which high- and low-momenum particles are 
indicated with thick and thin lines respectively. Explicitly the contribution of the diagram 
in Fig. 5 can be written as 

2TT5{kl -m'^)V{ki,k2\k[,k'2)V{ki, k2\k'l, 4')2Im$2(p; k[,k'2; k'(, k'^) (26) 

Here V is the relativistic (Bethe-Salpeter) potential and $2 corresponds to the lower blob 
with all intermediate states of A — 2 nucleons. 

The high-momentum part includes the first factor, potentials and the 6 -function coming 
from the real nucleon number 2. We can take the latter factors out of the integral at the 
point which corresponds to the initial nucleons at rest, that is at k'^ = k" = p, j = 1,2 and 
p = (m — e, 0). Then the high-momentum part separates as a factor, so that we get 

2Im^{p,ki) = 2Tr6{{ki - 2pf - m^)A^{A - l)w2H{p,ki). (27) 

Here the high-momentum factor H is 

Hij,M) = Hik!)^Y^i^^^!^^ (28) 

{H depends only on kf since the product pki can be related to it from the condition fef = 
{2p — fci)^ = m^). Factor W2 is just a constant which accumulates all the information from 
the nucleus: 

r ^ d'^k' k" 

A\A-1)W2 = j n ^^^^^\ 2T:)H\k'^ + k'2-kl-k'i)2lui^2(p;k[,k'2-,k'lk'i) (29) 

The calculation of W2-, is more or less standard. Similarly to Im$, the imaginary part of 
$2 can be presented in the form 

2lm^2{p-MM-Mk2) = j d'iiidVd^fe^^^+'''(*^i-^/2)-if(fci-if/2) 

{A\T{n\r + r/2)N\R - r /2)}T{N{f /2)N{-f /2)}\A) (30) 

where K = k\ + k2 = ki + k2 and, as before, N{N^) are the annihilation (creation) operators 
of the nucleon field and |^) is the nucleus state at rest. Putting this representation into (29) 
we find that the mometum integration puts all coordintes to zero: 

A\A - l)w2 = {A\{N^{0)fN\0)\A) (31) 
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To relate this expression with the standardly defined quantities, we note that all particles 
may be considered non-relativistic. So we may pass to non-relativistic states and operators. 

\A) = ^/2A^\AMR), N{N\r) = ^NMR{N],j,){r) (32) 

v2m 

Thus in terms of non-relativistic quantities we get (suppressing subindex NR) 

A\A-l)w, = ^{A\{N\^)fN\Q)\) = ^^^^^-^ j ndV,|V^(^)(0,0,r3,...rA)|^ (33) 

1=3 

where if)^^^ is the non-relativistic wave function of the target nucleus. It appears at zero 
distance between the 1st and 2nd nucleons, which is physically natural and corresponds to 
a correlation between them. To express (33) via better known quantities, one may take into 
account the translational invariance by presenting 

|V(^)(ri,r2,r3,...rA)P = j d?Tp^^\v^ + r,V2 + v, ...vj^+v) (34) 



with the nuclear p-matrix normalized as 

A 



(35) 



/ ndWHri,-rA) = l 
In terms of this /9-matrix we find from (33) 

W2 = ^Jd^rp(^\r,r) (36) 



where 



A 

P(ri,r2)= / l[d'np^^\ri,r2,rs,....rj^) (37) 
•' i=3 

and is normalized to unity after integration over ri and r2- The form (34) admits presenting 
p^^^ as a product of single-nucleon distributions. In this approximation 



^ 2m 



(iV(p(^)(r))2 (38) 
and thus is directly expressed via the nuclear density. 



4.2 Distributions in x and 

According to (11), to find the distributions in z and k'j_ we have to integrate the obtained 
function $(p, k) over This integration is trivial due to the 5-function in (27). It fixes the 
value of k'^ to be 

= zrr?^^ - kl^ - (39) 
2-z 2-z ^ ' 
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One finds that p{z, k\) is different from zero in the region (see Appendix 2) 

kl>A''(2-z-^^-m\z-lf (40) 

This condition is operative only if we take into acount the binding energy /m. Otherwise 
it is satisfied at all k_\_ and 1 < z <2. 

The internucleon relativitic potential can be taken as a sum of one-boson exchange con- 
tributions 

V{k,M\k[A) = Y.T?T^ (41) 

i f^i ^ 

where q = ki — k'ils the momentum transfer. In our case k'l = k2 = p, ki = k and ^2 = 2p — k. 
So one easily finds 

q^ = \{kl-m') + A^ (42) 
The distribution in z and k\ is obtained directly from (11):. 

^ = "^M^""^^'^ = "^ 4.(2'-.) [m^-k^f 5 (M? - g')(M? - g^) ^''^ 

where k"^ and q^ are expressed via z and by Eqs. (39) and (42) and W2 is given by (38). 

The remaining task is to integrate over k^ to obtain the final distribution in z, which is 
straightforward. Rather than discuss it we proceed to the realistic spinor case. 

4.3 Spinor nucleons 

A rigorous treatment of the spinor nucleon case requries knowledge both of the spinor struc- 
ture of the virtual photon-nucleon amplitude and the nucleus wave function, which is far 
beyond present (and future) possibilities. For this reason it is natural to use spin-averaged 

quantities for both quantities. More concretely we separate the integration over the virtual 
nucleon momentum and present the contribution from the diagram in Fig. 1 (right-hand side 
of Eq. (4)) in the form 

D = j d^A;Sp{(m + k)U{q, k){m + k)L{p, k)} (44) 

where U{p,k) is the upper blob (virtual 7-N amplitude) and L{p,k) is the lower blob (with 
factor l/{m? — fei)^ included). Our approximation is then to take 

Sp{(m + k)U{q, k)(m + k)L{p, k)} ~ ^Sp{(m + k)U{p, A;)}Sp{(m + k)L{p, k)} (45) 

which amounts to taking spin averaged values multiplied by the number of the spin states 
(two). After this averaging is made, all our formulas derived for the scalar case remain valid 
also for spinor nucleons. 
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Let us apply this recipe for the contribution from the two-nucleon correlation, correspond- 
ing to Fig. 5. In the spinor case instead of Eq. (26) we find for the function 2Im$(p, ki): 

2Im^p,h) = -^-^^^J Y[^-^^{27r)''d\k[ + k',-k'l-k'^)2^d{kl-m') 

Sp {2Im$2(p; k[, k'^-k'i. k'^){m+k'Mm+k2)2C{ki, k2\k[, k'^\k'i, ^^0(^+^1)1(^+^2)2} (46) 
where C is the two-nucleon correlation blob: 

C{k^, k2\k[, k'^lk'l, 40 = V{k[, k'^lki, k2){m + k){m + k2)V{ki, k2\k'l, 40 (47) 

and V is the relativistic interaction. Both C and V arc matrices in in spinor indeces of both 
nucleons. Operators (m + /c)i(2) act on the spinor indeces of 1st (2nd) miclcon. The trace is 
to be taken over spinor indeces of both nucleons. In deriving (46) we have applied (45) to 
average over the spin states of the active nucleon (momentum ki). Note that the factor 1/2 
is to be included in the definition of the spin- averaged nucleon structure function F in Eq. 
(5). 

Now we use the fact that the momenta k[,k2, k'{ and k2 are non-relativistic and also apply 
the recipe (45) two more times to sum over spin states of the nucleons 1 and 2 before the 
left interaction and after the right interaction in Fig. 5 to obtain a representation identical 
to (27) where now 

H{kj) = ^Si>{{m + p)iim + p)2C{ki,2p - ki\p,p\p,p)}, (48) 

and 

A^A - i)w2 = / n ^^^{2^mk[ + 4 - k'l - 40 

Sp { 2 Im $2 {p; k[ , 4 ; 4', 40 (m + p) 1 (m + p)2 } (49) 
In the latter expression in the non-relativistic approximation we may write 

Sp { 2 Im ^2 {p; k[ , 4 ; k'l, k'^){m + p) 1 (m + p)2 } 
~ ^Sp {(m + 4)(m + 4) 2Im$2(p; 4,4; 4', k'^){m + k'Dim + 4)} 

Sp {A\T{N^R + r/2)N^R - r/2)}T{N{f/2)N{-r/2)}\A) (50) 

where now iV(iVt) are 4-spinor operators of the nucleon field. In the non-relativistic limit their 
first two components reduce to non-relativistic spinor fields for the annihilation (creation) of 
a nucleon with a given spin component at a given point, the second two reducing to zero. As 
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in (30), the integration over the momenta in (50) puts all coordinates to zero and we obtain 
similarly to (33) 

(T,(T2 = 1 

= J n^'^^n E KU,.....(0,0,r3,...rA)P (51) 

i=3 j=laj=l 

where cr's denote non-relativistic spin projections of the nucleons. Obviously the integral on 
the right-hand side has the same meaning as in the scalar case: the total probability to have 
the two active nucleons at the origin irrespective of their spins. Presenting the spin-summed 
square modulus of the nuclear wave function in the form (34) one arrives at the same formulas 
(36) or (38) as for the scalar case. 

So in the end the only difference from the scalar case is the expression (48) for the function 
if as a spin average over the two-nucleon correlation contribution given by Eq. (47) in terms 
of the relativistic internucleon interaction. 



5 Numerical results 

Calculation of the cumulative nuclear structure function clearly splits into three steps. First 
the density p{x, k) has to be determined from the given relativistic internucleon potential V. 
Second this density has to be intcrgratcd over the transverse momenta to find the resulting 
distribution in x. Finally one should convolute this distribution with the known structure 
function of the nuclcon. 

To determine p(a;, k) one first should calculate the two-nucleon correlation factor H ac- 
cording to Eqs. (47) and (48). We borrow the relativistic internucleon potential V from [5], 
where it is represented by a sum of vr, rj, a, oj and p exhanges, modified by form-factors to 
cut the high momentum contribution and adapted to off-shell nucleons. In fact four different 
parametrizations of the potential are presented in [5] . They are given in Appendix 1 together 
with the concrete form of various contributions. As we shall see different parametrizations 
give practically the same results for the cumulative densities and structure functions. 

The calculation of H is rather tedious due to many different terms and relativistic traces. 
Some intermediate formulas are given in Appendix 2. 

The nuclear factor was calculated from the standard Woods-Saxon nuclear density. 
Obviously it is of the order l/A and has a dimension of mass squared. Correspondingly in 
Fig. 6 we show its values multiplied by {A — l)/w? for different nuclei. The resulting values 
steadily grow with A. Their ^l-dependence does not correspond to A°^: one finds that in this 
parametrization a goes down from 0.3 to 0.15 as A rises from 20 to 200. 

Intergrating the found p(x, k) over the transverse momenta we found the distribution 
p{x). In Fig. 7 we show its values for Carbon with two potential models from [5]: lA and 
IIA. One observes that up to x ~ 1.3 both potentials gives the same p{x). At higher x model 
lA gives values for p{x) higher by factor 1.5 ^ 3, which is quite a small difference on the 
cumulative scale. 
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The resulting cumulative structure function of Carbon are presented in Figs. 8 and 9 at 
= 3,30 and 1000 (Gev/c)^ for potental models lA and IIA respectively. For the nucleon 
structure function F{x,Q'^) the GRV LO paramctrization has been taken from [6]. The 
results are quite similar for the two potential models, both in form and magnitude. One 
observes that the cumulative structure function turns out to be rather weakly dependent on 
Q^: at moderate x it falls with by a factor of the order 4, as rises from 3 to 1000 
(GeV/c)^. As to the x dependence it is evidently non-exponential, although universal in Q^, 
the slope of F^^\x) changing from 10 to 40 as x rises from 1.2 to 1.7. 

Comparison with the experimental data from [2] is presented in Table 1. One observes that 
at a; = 1.05 and 1.15 our results lie above the data by roughly 100% and 50% respectively. 
This disagreement is not so dramatic considering the steep fall of the structure function 
beyond x = 1. One has also take into account that these values of x are too close to unity 
for our formalism to work well. 

Table 1. Structure Function of Carbon (x 10^) 
in different potential models of [5] 



X 


(GeY/cf 


lA 


IB 


Ila 


IIB 


Exp 


1.05 


61. 


162. 


153. 


126. 


153. 


64.1±7.4 




85. 


150. 


142. 


116. 


142. 


47.8±5.2 




150. 


133. 


126. 


103. 


125. 


45.8±5.1 


1.15 


61. 


39.5 


40.7 


30.9 


39.7 


<27.0 




85. 


36.4 


37.6 


28.5 


36.7 


<15.1 




150. 


32.0 


33.1 


25.0 


32.3 


<22.3 


1.30 


61. 


3.65 


4.15 


2.71 


3.89 


<6.9 




85. 


3.35 


3.81 


2.48 


3.56 


<5.4 




150. 


2.90 


3.30 


2.15 


3.09 


<3.8 



6 Conclusions 

We have developed a method which allows to calculate the nuclear structure function at 
a; > 1 in terms of nucleonic degrees of freedom with full relativistic kinematics. The method 
avoids solving the complete relativistic many-body problem and relates the structure function 
directly to the high-momentum asymptotics of the relativistic potentials. While for the 2- 
body correlations it can be considered as a simplification of the solution of the Bethe-Salpeter 
or Gross equations, for the correlations involving more nucleons it seems to be the only 
realistic approach. 

The results of our calculations confirm that the nuclear structure function at x > 1 
depends on only weakly. Its x-behaviour cannot be well desribed by an exponential. 
Comparison to (scarce) experimental data for C at x = 1.05 and 1.15 shows that our results 
have the same order of magnitude, lying above the experimental data by 100% and 50% 
respectively. This disagreement may be completely explained by the comparatively low values 
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of X where our formalism is not supposed to work well. One should wait for more data at 
higher x to draw final conclusions as to the validity of the chosen picture and the potentials. 
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9 Figure captions 

Fig. 1. The impulse approximation for the nuclear structure function. 
Fig. 2. The impulse approximation for the nuclear form-factor. 

Fig. 3. The asymptotics of the Bethe-Salpeter function from an iteration of the acting po- 
tential. 

Fig. 4. The two-nucleon (a) and three-nucleon (b) correlations in the nucleus. 
Fig. 5. The contribution from the two-nucleon correlation (Eqs. (26) and (46)). 
Fig. 6. The nuclear factor (^4 — l)w2/'m'^ as a function of A . 

Fig. 7. The distribution p{x) for ^^C in the cumulative region for potential models lA (upper 
curve) and IIA. 

Fig. 8. The cumulative structure function of ^^C with the potetial model lA. Curves from 
top to bottom correspond to = 3, 10,30, 100 and 1000 (GeV/c)^. 
Fig. 9. Same as Fig. 8 with a potential model IIA. 
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